Optimization of transport protocols with path-length constraints in complex networks 
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We propose a protocol optimization technique that is applicable to both weighted or unweighted 
graphs. Our aim is to explore by how much a small variation around the Shortest Path or Optimal 
Path protocols can enhance protocol performance. Such an optimization strategy can be necessary 
because even though some protocols can achieve very high traffic tolerance levels, this is commonly 
done by enlarging the path-lengths, which may jeopardize scalability. We use ideas borrowed from 
Extremal Optimization to guide our algorithm, which proves to be an effective technique. Our 
method exploits the degeneracy of the paths or their close-weight alternatives, which significantly 
improves the scalability of the protocols in comparison to Shortest Paths or Optimal Paths protocols, 
keeping at the same time almost intact the length or weight of the paths. This characteristic ensures 
that the optimized routing protocols are composed of paths that are quick to traverse, avoiding 
negative effects in data communication due to path-length increases that can become specially 
relevant when information losses are present. 

PACS numbers: 89.75.Hc, 89.20.Hh, 05.70.Ln 



I. INTRODUCTION 

Communication systems are an essential element of the 
modern world, and represent some of the most notice- 
able settings of transport problems in complex structures. 
The Internet, telephone systems, power grids, etc., pro- 
vide ways to transport information, signals or resources 
from one part of the globe to another. Such commu- 
nication is done by means of highly structured systems 
that can be modeled by networks, i.e., collections of el- 
ements corresponding to nodes, links that interconnect 
them, and weights on the links symbolizing connection 
capabilities [IH5]. 

In many examples of communication systems, discrete 
subsets of information (packets) are sent from a source 
to the appropriate destination. Typically, a multitude of 
users simultaneously perform similar tasks between many 
different source-destinations pairs, and thus the network 
at any given time is populated by numerous traveling 
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data packets. In order for packets to find their way, rout- 
ing protocols are put in place which determine how the 
packets travel from source to destination. Optimization 
of the performance in communication systems can thus 
be achieved in two non-exclusive ways: improving the 
structure of transport networks [B] , and finding optimal 
routing protocols [7] . Both approaches have been the fo- 
cus of a great deal of research on this subject in recent 
years [6lI25j. A number of important results have been 
obtained, elucidating the usefulness of so-called scale-free 
networks (7| [T^] , the vulnerability of networks to random 
or targeted failures of some of its elements [H] , and the 
relevance of path-lengths in the efficiency of communica- 
tion [TSJ. 

When the overall traffic in the network is very large, 
the system may get compromised in its capability to de- 
liver information. The origin of this problem could be 
twofold: either the connections of the network are unable 
to tolerate such a high traffic flow or, more commonly, 
some components (nodes) get overloaded, i.e., congested. 
Congestion occurs when the nodes become heavily used 
acting as bottlenecks for the network communication. By 
establishing clever routing protocols, this congestion can 
be diminished or, for moderate rates of overall traffic in 
the system, eliminated. The assessment of node work- 
load requires thus the introduction of a new quantity. 
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Betweenness centrality (betweenness in short) has been 
proposed as such a measure [BJ [17] , as an estimate of how 
much use a given node can expect to receive. Specifically, 
for a set of paths in a network, the betweenness of a node 
is defined as the number of paths visiting that node. 
Therefore, a given routing protocol induces a between- 
ness on the nodes of the network. It has been pointed out 
that nodes with large betweenness are the first to form 
packet queues when traffic is large enough [BJ. If this 
queue is not eliminated, then the travel times of pack- 
ets from source to destination increase without bound, 
hence, rendering the node (and potentially the network) 
congested. 

On the other hand, packet delivery is an imperfect pro- 
cess since it can be subjected to packet loss or corruption. 
When congestion occurs, even entire queues of packets 
could be eliminated while waiting at overloaded nodes. 
Thus, when attempting to optimize a routing protocol, 
care must be taken not to increase the path-length to a 
level that would elevate the likelihood of packet loss be- 
yond critical limits. Such considerations have led to the 
formulation of models such as Limited Path Percolation, 
which exhibit a phase transition when the length of paths 
increases beyond a specific tolerance level [T5] . 

In this article, we focus on the problem of optimizing 
the transport rules for scale-free networks with proper- 
ties similar to those of computer networks. However, 
in contrast to a number of other studies, we optimize 
the protocols to reduce the likelihood of congestion but 
keeping path-lengths as close as possible to the short- 
est, thus guaranteeing efficient communication paths. We 
find that by applying ideas borrowed from Extremal Op- 
timization (EO) |26| . which uses local information to ap- 
proach the optimal state of a system, we can achieve both 
goals of congestion avoidance and short path- lengths. 

The article is structured as follows: Sec. |H| explains the 
philosophy behind our optimization heuristic, the algo- 
rithm to find optimal routing protocols, and our network 



model. Sec. Ill quantifies betweenness of nodes as a con- 
sequence of the optimized routing protocol. Sec. |IV] fo- 
cuses on the path-lengths created by the new, EO-based, 
routing protocol, and the advantages that it affords re- 
garding maintaining the path-lengths within a constant 
factor of the original paths. Sec. [V] extends our methods 
to weighted networks, and finally, Sec.|Vl]is dedicated to 
the conclusions. 



II. OPTIMIZATION METHOD 

To define the betweenness precisely, we need first to 
specify a routing protocol, which is a set of paths C — 
{pij} between all pairs of nodes i and j of the network. 
We define the length of each path lij as the minimum 
number of links needed to be crossed when going from i 
to j following the path pij . We adhere to the convention 
of considering only one path for each pair of nodes [20] , 
as the length degeneracy of the paths will be exploited 



by the method to optimize the protocol. The rules by 
which these paths are chosen can be numerous. For in- 
stance, one typical choice for C in the literature is the set 
of shortest paths between nodes (SP). The betweenness 
bi of a node i is then defined as the number of paths of 
C that visit node i, without counting those paths that 
begin or end in i, (This choice is somewhat arbitrary, 
and other choices are also possible.) Note that in the 
literature the word "betweenness" sometimes is associ- 
ated exclusively with the number of paths crossing a node 
within the framework of the SP protocol. We are using 
this concept in a more general context since it can re- 
fer to any routing protocol, which has to be given. In 
fact, betweenness values depend on both, the structure 
of the network and C For a given network, changing the 
routing protocol typically changes bi over the network 
except in very special cases such as when there is only 
one possible protocol choice (e.g., a star network). 

Given a network and a protocol C, it has been shown 
that the first node to begin accumulating a queue, say 
m, satisfies that b m > bi for all nodes i ^ m in the net- 
work [BJ. To explain this point, consider a network of size 
N in which data packets are produced in each node at 
random with a rate < 7 < 1. In the simplest version 
of the problem, the destination of each packet is random 
too, and thus, each node can produce a packet to any of 
the other N—l nodes of the network with equal probabil- 
ity. In this case, it was found that the number of packets 
in node i is, on average, r yb i /{N — 1) [BJ. Without loss 
of generality, we assume that each node can processes 
(route through) a single packet per unit time, and thus, 
the first node to congest is the one with the maximum 
betweenness value, & max . We can also rephrase this point 
by indicating that the critical rate of packet production 
would be 7 C = (N— l)/b max before some node of the net- 
work congests [BJ. This relation shows the relevance of 
designing routing protocols with a £> max as low as possi- 
ble in order to prevent network congestion for low 7 val- 
ues. The lower 6 max is, the higher the number of packets 
can be (higher 7) before the network congests. There 
is, however, a price to pay in this strategy because the 
reduction of 6 max is usually attained at the expense of 
enlarging the path-lengths. Under some circumstances, 
such as when packets can get lost, enlarging significantly 
the paths may not be the best option. The model that 
we propose next aims at bridging the gap between these 
two opposing forces: we explore how much 6 max can be 
reduced while keeping the paths as close to their shortest 
version as possible. 

We generate the networks for the analysis with the con- 
figuration model [371 [21] that allows us to build graphs 
with arbitrary degree sequence but free from degree- 
degree correlations. The degree of a node i corresponds 
to its number of connections and is usually represented 
as hi. In this work, we will focus on power-law de- 
gree distributed networks, P(k) ~ k ~ A , with exponent 
A = 2.5 to facilitate the comparison with previous stud- 
ies [2BJ dU [23] . Nevertheless, we have also performed the 
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FIG. 1: (Color Online) The sketches (a) and (b) illustrate the 
way in which the method optimizes the protocol by selecting 
another path that must be either of the same length or the 
shortest alternative to the original one. In panel (c), the prob- 
ability of finding alternative length-degenerate paths, P n >i, 
for the Shortest-Path protocol is displayed as a function of 
the network size N. In the inset: histogram for the number 
n of degenerate alternatives for each path in the SP proto- 
col. The different curves correspond to the following network 
sizes: N — 100 (circles), 500 (squares), 2000 (diamonds) and 
10000 (triangles). 



stressing that, the constraint on the length of the alter- 
native paths in practice means that if another choice of 
equal length to the original path exists, the path does not 
increase in length. The method takes advantage of the 
possible degeneracy of the paths with respect to length 
to change the protocol and ease the stress on the node(s) 
with the largest traffic. Note also that the method keeps 
no long-term memory, the search for alternative paths 
due to sequential expulsions of a path from different 
nodes do not accumulate. The algorithm generates thus 
protocols with paths limited to SP configurations or the 
shortest alternatives to them avoiding a single high used 
node. 

The potential utility of exploiting the SP degeneracy 
and its importance in our networks can be seen in Fig- 
ure [I] We have plotted in the panel (c) the probability 
of finding a degenerate alternative to the shortest path 
between any pair of nodes i and j, P n >i{N), as a func- 
tion of the network size N. P n> ±(N) consistently grows 
with TV, becoming higher than 65% for the largest graphs 
that we have considered N = 10000. In the inset, the 
distribution of the number of different degenerate paths, 
P(ri), is also shown for four network sizes. Note that 

The philosophy of our method is in the spirit of Ex- 
tremal Optimization (EO) [25], attacking the worst ele- 
ment of the system in an attempt to improve the global 
performance. In the following, we will study how its ap- 
plication affects the different features of the protocols. 



III. BETWEENNESS DISTRIBUTION 



analysis on networks with A = 4.5 finding no qualitative 
difference in the final results. 

The optimization algorithm starts from the SP proto- 
col: Once the network is built, the first step is to calculate 
the shortest paths between all nodes of the network using 
Dijkstra's algorithm |29j . Then, in an iterative way for 
as long as desired, we repeat the following steps: 

• The node with the highest betweenness m is se- 
lected. 

• One of the paths, pij, passing through m is chosen 
at random; we refer to its initial and final nodes as 
i and j (with i =/= m =/= j). 

• A path between i and j is searched. The new path 
must be the shortest alternative to not passing 
through m. In case that there are more than one 
alternative paths degenerate in length, one of them 
is selected at random. On the other hand, if there 
is no alternative to py, the path remains without 
change. 

A schematic illustration of these rules is shown in pan- 
els (a) and (b) of Figure I] In the sketch, a path pass- 
ing through the highest betweenness node, the square, 
is re-routed thereby reducing betweenness. It is worth 



We will discuss first the results concerning the be- 
tweenness of the nodes and how its values adapt after the 
optimization method is employed. After some iterations 
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FIG. 2: (Color Online) Betweenness distribution for a N = 
500 system. The (black) circles correspond to the stationary 
distribution reached with our optimization method while the 
(red) squares are for the initial SP protocol. In the inset, the 
stationary betweenness distributions of the method for other 
two network sizes: N = 350 and N = 700. 



4 



10 



10 r 



S 3 

£ 10 



10" 



10 



Protoc. 


a 


SP 


1.62 


HA 


1.48 


DAN 


1.18 


ER 


1.31 


Our 


1.3 


i i i i i 



20 



100 



N 



1000 



FIG. 3: (Color Online) The average maximum betweenness, 
bmax, as a function of the network size N for the optimized 
protocols. Lines are shown as a reference scaling with ex- 
ponent 1.3 (solid), 1.18 (dashed), and 1.62 (dotted). In 
the table, scaling exponents a obtained with different op- 
timization methods are listed. The initials correspond to 
Shortest-Paths (SP) [31], Hub- Avoidance (HA) [20], Danila's 
method (DAN) [ST], Efficient Routing (ER) [S3] , and our own 
method. 



of the method [3D] the betweenness distribution, P(b), 
reaches a stationary form that can be seen in Figure [2] 
Before describing it, it is important to remember that the 
SP protocol produces a power-law distribution for P(b) 
whose exponent depends on the exponent A of the de- 
gree distribution [TT] [TB] (see the curve of (red) squares 
in Figure [2]). The protocols generated by the optimiza- 
tion method show that the tail of P(b) collapses into a 
peak at high values of the betweenness. Also the area of 
low b suffers a slight variation, although the functional 
form of P(b) at intermediate b seems largely unaffected. 
A change in the size of the network displaces the position 
of the peak but not the quality of the effects observed in 
the main plot of Figure [2] 

The peak induced in P{b) by the optimization method 
at high values of b correlates with a decline of the maxi- 
mum betweenness in the network. As mentioned before, 
this betweenness, b max , is of crucial importance for the 
transport capability of a protocol, since it imposes an 
upper cutoff to the traffic that the network is able to 
sustain before congesting [DJ [T3J. The scaling of b max 
with the network size, N, can thus be seen as a mea- 
sure of the scalability of a routing protocol: the faster 
bmax grows with N, the less useful the protocol is for 
large networks. In Figure [3j the average b max in the 
stationary state of the protocols obtained with our op- 
timization method is depicted as a function of the num- 
ber of nodes in the network. We find a power-law in- 
crease with a functional form b max ~ N a and a our ps 1.3. 
This value of a must be compared with the results en- 
countered with other protocols or other protocol opti- 
mization methods. In the same Figure, we have also in- 
cluded a table with the exponents a for other methods. 




FIG. 4: (Color Online) In panel (a), path- length distribution 
for a network of size N = 1000. In the figure, the results of 
the protocols obtained with our optimization method can be 
compared with the initial SP. In panel (b), the ratio between 
the average length of the paths in the two protocols (ours and 
SP) is visualized as a function of the network size. 



It is specially interesting to compare with the original 
Shortest-Path protocol a SP ps 1.62 [3T] or with the best 
of the optimization methods listed (Danila's method) 
with a DAN « 1.18 [1UI12]. Our optimization method, 
without deviating substantially from the SP protocol in 
the length of the paths, produces an acceptable value of 
the exponent a. 



IV. PATH-LENGTHS 

The near-invariance of the path-length between SP and 
our optimized protocol, and the path-length distribution 
itself, deserves further attention. Here we will focus on 
studying how the length of the paths in the protocol mu- 
tates when the optimization method is applied to the 
protocol. In Figure [4j we have plotted the path-length 
distribution P(£) for the original Shortest Path proto- 
col (red squares) and for the stationary regime of the 
protocols obtained with our optimization method on a 
network of size N — 1000. As can be seen, the distri- 
bution is slightly displaced to larger values of I but its 
shape remains essentially invariant. In order to quantify 
the relevance of the change of P{i) with the system size, 
we have represented in panel (b) of Figure [4] the ratio 
between the average length of the paths in the protocols 
obtained with our method and that calculated with the 
SP. The curve is monotonically decreasing with N, mov- 
ing very slowly towards unity (or a value just above). 

Similar to b max , limiting the scaling of the average 
length of the paths, (£), with the system size is also an 
important design feature. A protocol that elongates the 
paths unnecessarily may not be efficient for network com- 
munications. Generally, there exists a nonzero probabil- 
ity of losing a packet in every communication between 
two servers, and the longer the paths become, the higher 
the probability of missing information. Furthermore, if 
data loss reaches a point in which most of the paths are 
not functioning, the network suffers a disruption similar 
to a percolation transition. In fact, this phenomenon has 
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been studied recently and it is known as Limited Path 
Percolation [15]. In Figure [5j we show how {£) in the pro- 
tocols obtained using our method behaves with increasing 
network sizes. In scale-free networks with A = 2.5, the SP 
protocol (£sp) is expected to grow as log(log(7V)). Wc 
cannot numerically test this formula for many decades 
in N for our optimized protocols, but within the limited 
range of values that we can explore the relation holds. 
Such scaling is important because, as mentioned above, 
the optimization of the protocols can in some cases im- 
ply a trade-off between reductions in b max and escalating 
path-lengths. For instance, Danila's method [3T], which 
exhibits the best scaling for b max versus N of the methods 
shown in Figure [3j leads to a (£) ~ log(iV) [22], consid- 
erably faster than the log(log(iV))-scaling of our method 
(or of the SP protocol). This difference in the scaling 
can be relevant if the packet communication has losses in 
long paths (or in paths longer than the SP). 

In general, consider a protocol that is tuned to scale 
with the typical structural path lengths of the network, 
which, to be concrete, we associate with the most prob- 
able length £* SP of the length distribution P SP (l). We 
imagine that the real system uses a protocol C when 
uncongested which is, for our purposes, the best practi- 
cal protocol for the network of interest (even if it is not 
strictly optimal) . For a power-law network with exponent 
A = 2.5, t = D t sp ~ D loglog(iV), where D Q > 1 is 
a numerical constant. Now, when higher rates of packet 
creation appear on the network, to avoid congestion we 
consider introducing two new protocols £5 and Cl, the 
first one with typical path lengths l* s ~ 8si* {5s > 1 
another constant), and the second with l* L = F(t f ) 1 with 
F(x) > x 1+e , e > 0. In the thermodynamic limit, with 
mean field approximation, r, the ratio of new typical path 
length to old typical path length in £5 remains constant, 
as we can sec from 
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FIG. 5: Average path- length, (l), of the protocols we obtain as 
a function of the network size. It can be seen that a log(log())- 
scale in panel (b) provides a better fit to our simulation results 
than a pure logarithmic scale in (a). 



where the last relation occurs for N — > 00. In Ref. |15j . 
the limited path percolation transition was shown to be 
intimately related to r. For our purposes, it means that 
if a protocol keeps r finite as N — > 00, occurring in the t$ 
case, there is the possibility to have a protocol that can 
be tuned to maintain at least a level of communication 
between node pairs. Note that the C$ protocol can be 
written as F{x) = x. On the other hand, if r diverges, 
as in tl, as the system size grows, the typical pair of 
nodes becomes increasingly distant in the relative sense, 
and thus cannot in general remain in communication, or 
in a more formal sense, protocols such as Cl may not be 
scalable [55] , 

The protocol we propose here has the behavior of Cs 
and thus, as in limited path percolation, can have a per- 
colation threshold to a connected regime. Furthermore, 
the fact that we study scale-free degree distributed net- 
works with small A indicates that from the standpoint 
of limited path percolation, there is always a connected 
regime. In addition, when 5 is close to 1, that the pro- 
tocol is very robust as most pairs of nodes remain con- 
nected after the introduction of the optimization of b max . 
These conclusions are also supported by results displayed 
in Figures [4j [5] and [6j 

Furthermore, as explained above, a path obtained with 
our method can be only in two states: either the SP 
configuration or, if they were expelled from a high used 
node and no lcngth-dcgcncrate alternative was available, 
in the shortest configuration that does not pass through 
that node. Such two state configuration recalls physical 
systems in which the components can be in the energetic 
ground state or being activated by thermal (stochastic) 
fluctuations to higher energy levels. To confirm whether 
this parallelism is valid, we have plotted in Figure [6] the 
cumulative distribution (C > (R) — dR' P(R')) of the 
ratios R — l s tat/£sp, where £ sta t corresponds to the 
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FIG. 6: Color Online) Cumulative distribution of the ratio 
between the length of the paths obtained with our protocol 
optimization method and those of the corresponding original 
SP. The distribution becomes step-like due to the discreteness 
of the length. 
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length of each path in an stationary instance of the pro- 
tocols obtained with our optimization technique and £sp 
is the corresponding length of the shortest path. As can 
be seen, the probability for the paths to be outside the 
ground state (SP), in an "excited state", decays as an 
exponential function with a characteristic ratio that de- 
pends on the graph (whether the graph is scale-free or 
not, and, if it is, also on the exponent of the degree dis- 
tribution). In the example of Figure |6j the characteristic 
ratio is around 1.55, which makes it extremely unlikely 
for any path to suffer a large stretching out of its SP 
length. This exponential decaying C>(i?) resembles the 
Boltzmann factor in which R — 1 plays the role of the 
energy. 



V. WEIGHTED NETWORKS 

The introduction of disorder in the network connec- 
tions deeply modifies the previous scenario. However, 
connections of divers quality can be expected in may real- 
world situations. A classical example within the context 
of transport is the Internet and the variability on band- 
width of the connections between servers. Mathemati- 
cally speaking, the quality of a connection can be repre- 
sented by a scalar attached to each link that is known as 
its weight w. The presence of weighted links alters the 
definition of shortest path 32 34J . Instead of "short" in 
the sense of number of links between origin and desti- 
nation, it may be important to find the path with the 
lowest weight along the way (or the largest, depending 
on the quality that the weight conveys) . These paths are 
often referred to as optimal paths in the literature (see, 
for instance, [55]). 

Protocols based on Optimal Paths (OP) also have a 
very particular scaling with the network size [13j ■ De- 
pending on the level of randomness on the weights of 
the connections, the system can fall into a strong- or 
a weak-disorder regime. The difference between these 
two regimes is that in the strong-disorder case the fluc- 
tuations on the weight are so large that the accumu- 
lated weight along each paths is controlled by the edge 
with the highest weight, while in weak-disorder the "re- 
sponsibility" for the path's overall weight is distributed 
among all the links along the path [33 [37] . In the two 
regimes, the scaling of b max , and of the average weight 
of the paths (w pat h}, with N changes with respect to 
unweighted graphs [2SJ. Here we will focus only on the 
weak-disorder regime since it is the only one in which 
it makes sense to search for an alternative protocol to 
the OP. The cost of doing so in the strong-disorder limit 
would diverge with the size of the graph. 

One of the positive aspects of our protocol optimiza- 
tion is that its generalization to weighted graphs is 
straightforward: at each step, a path passing through 
the node with the highest betweenness is chosen to be 
redirected. The path is then exchanged by its alterna- 
tive of lowest cost that does not pass through that node. 
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FIG. 7: (Color Online) Dependence of b max with the system 
size for weighted graphs. The (red) circles are the results 
for the protocols of our method and the (blue) squares are for 
the original OP protocol. The weight distribution in panel (a) 
follows an exponential P(w) ~ e~ w ^ 2 , and in (b) the power- 
law P(w) ~ w~ 2 ' 5 , as displayed in the insets for networks 
of size N = 500. The continuous curves in both panels are 
a guide to the eye showing a visual comparison between the 
exponent values obtained by fitting the simulation results. 



As before, if there is no alternative, the path remains in- 
variant. In the following, we will study the performance 
of the stationary protocols produced by this method on 
Molloy-Reed graphs with A = 2.5 and with two possible 
functional forms for the weight distributions: either an 
exponential P(w) ~ e -w/w c or a power-law P(w) ~ w^ 13 . 
In order to keep the results comparable for both distri- 
butions, we fix the parameter w c to obtain a similar av- 
erage weight as in the power-law distribution. For the 
data shown, this will be w c — 2 and j3 = 2.5. 

The scaling of b max as a function of the network size 
is shown in Figure [7j for networks with both type of 
weight distributions. The application of the optimization 
method produces a significant improvement in b max with 
respect to the Optimal Paths protocol. For the range of 
system sizes of the figure, this entails an improvement 
by a factor 3 or even higher. Additionally, the exponent 
a of the optimized protocols is close to 1.4 while that 
of the original OP protocols was about 1.55. One can 
wonder about the price to pay for such an improvement. 
In networks with continuous values of the weights, any 
degeneracy in the weight of the alternative paths is ex- 
tremely unlikely. Strictly speaking, then, there is always 
a price to pay, although it could be minute. To answer 
this question, we have plotted in Figure [§] the average 
weight of the paths scaled with N for protocols obtained 
with our method and for the corresponding OP on both 
types of networks and both weight distributions. Sim- 
ilarly to the unweighted case, we find that the average 
weight is slightly higher but not extremely so. Also the 
ratio between the average weight of the paths of our op- 
timized protocols and of OP decreases for increasing net- 
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FIG. 8: (Color Online) Comparison between the average 
weight of the paths, (w pa th), in the original OP protocol and 
in the optimized protocols with our method. The curves in (a) 
correspond to weighted graphs with an exponentially decay- 
ing weight distribution, while those in (b) have a power-law 
weight distribution. 



work sizes. Furthermore, even though it is not shown in 
the figure, neither the shape of the distribution nor the 
weight of the paths is hardly altered by the optimization 
of the protocols. 



VI. DISCUSSION AND CONCLUSION 

We study a routing protocol optimization method that 
takes advantage of the degeneracy of the shortest paths 
to reduce the congestion in the nodes of the network. 
The paths between node pairs in the optimized proto- 
cols can only be in two states: they can be either one 
of the shortest-paths or the closest alternative avoiding 
highly used nodes. We use ideas coming from Extremal 
Optimization to guide our algorithm, which proves to be 
an effective technique. The resulting routing protocols 
show an appreciable reduction of the maximum node be- 
tweenness, 6 max , while keeping the path-length distribu- 
tion similar to the lower bound provided by the SP pro- 
tocol. For the networks that we have analyzed 6 max in- 
creases with N as a power-law with an exponent a ~ 1.3. 
That is only slightly higher than the best values reported 
in the literature, and, in contrast to those methods, the 
average path-length appears to retain the log (log (N))- 
scaling of the SP protocol. We show that methods which 
do not conserve a constant ratio of path-length with the 
topological shortest paths end up in a disconnected phase 
of percolation with limited paths, which is a possible in- 
dicator that on networks with packet loss such methods 
may not be scalable. 

Our method is easily generalized for weighted net- 



works. We have studied two examples of weighted net- 
works in the weak-disorder regime and found significant 
reductions in 6 max with respect to the Optimal Paths 
protocol as well as in the scaling exponent with which 
& max grows with N. 

As we implemented the optimization method-a pro- 
cedure in which we keep record of all the paths passing 
through each node-, the numerical performance is not 
optimal. Indeed, we were able to study a range of net- 
work sizes wide enough to characterize the scaling of the 
optimized protocols properties, but not as large as to al- 
low for large network analysis. Our aim in this work was 
to explore the behavior of the optimization method and 
to ascertain whether the resulting routing protocols have 
interesting features concerning the size scaling and be- 
havior of &max and the path-length distributions. Indeed, 
the simple fact that we have shown that it is feasible to 
perform such dual optimization is an important step for- 
ward. Furthermore, we have found that this two param- 
eter optimization is important because without taking 
into account both b ma x and path-length, the protocols 
that simply reduce b max may be less scalable. 

Therefore, the possibility of searching for efficient nu- 
merical methods to implement the optimization algo- 
rithm remains open for future work. A possible solu- 
tion, depending on the structure of the network, can be 
the implementation of a local search technique for the 
alternative paths. This local method could be highly ef- 
ficient if the clustering (the density of triangles) in the 
graph is high as typically occurs in geographically em- 
bedded or social driven networks. Locally tree-like struc- 
tured networks, such as those considered here, typically 
require global rearrangements between paths of nearly- 
degenerate length but hardly any overlap, as illustrated 
by Fig. la-b. Maintaining the paths of the routing pro- 
tocol short can be crucial in situations such as when the 
system suffers packet losses, a circumstance that seems 
to be relatively frequent in practice. Our optimization 
method represents in such cases a useful tool that offers 
a balance between betweenness reduction (lower conges- 
tion) and restricted path-lengths (reduced losses) in rout- 
ing protocols. 
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